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We study the effects of rapid temperature and magnetic field changes on the non-equilibrium 
relaxation dynamics of magnetic vortex lines in disordered type-II superconductors by employing an 
elastic line model and performing Langevin molecular dynamics simulations. In a previously equili¬ 
brated system, either the temperature is suddenly changed, or the magnetic field is instantaneously 
altered which is reflected in adding or removing flux lines to or from the system. The subsequent 
aging properties are investigated in samples with either randomly distributed point-like or extended 
columnar defects, which allows to distinguish the complex relaxation features that result from either 
type of pinning centers. One-time observables such as the radius of gyration and the fraction of 
pinned line elements are employed to characterize steady-state properties, and two-time correlation 
functions such as the vortex line height autocorrelations and their mean-square displacement are 
analyzed to study the non-linear stochastic relaxation dynamics in the aging regime. 
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I. INTRODUCTION 

Type-II superconductors have been the focus of many 
theoretical and experimental investigations since their 
discovery. From the fundamental physics perspective, 
these materials represent a class of systems with many 
competing energy scales and thus a remarkably rich vari¬ 
ety of thermodynamic phases, as well as intriguingly com¬ 
plex transport properties [l|. Applied science and engi¬ 
neering research on type-II superconductors has focused 
primarily on their optimization for technological appli¬ 
cations, especially in magnetic fields. The most promi¬ 
nent desirable property of superconductors is of course 
the absence of current flow dissipation in these materi¬ 
als. However in the mixed phase of type-II supercon¬ 
ductors, magnetic flux enters the sample in the form of 
quantized vortex lines, and external driving currents gen¬ 
erate a Lorentz force that in turn causes these vortices to 
move, creating Ohmic dissipation in the system. To pre¬ 
vent flux flow and thus restore dissipation-free transport, 
effective pinning mechanisms of flux lines to material de¬ 
fects are required [I|. 

Different types of defects are utilized as pinning cen¬ 
ters for magnetic vortices, predominantly uncorrelated 
point-like disorder and extended linear (columnar) de¬ 
fects. Point defects are naturally occurring, e.g., in ce¬ 
ramic high-Tc materials in the form of oxygen vacancies, 
but can also be artificially introduced, for instance by 
electron irradiation Q . The presence of weak point pin¬ 
ning centers destroys the long-range crystalline order of 
the low-temperature Abrikosov flux line lattice in the 
disorder-free systemMo form either a genuine disordered 
vortex glass phase [J-Q or a Bragg glass state that is 
characterized by quasi long-range positional order 
[T5| . The thermally induced first-order melting transition 
of the vortex lattice at elevated temperatures (l^-[l^ is 


thereby replaced by a disorder-driven continuous phase 
transition between frustrated (‘glassy’) low-temperature 
states and a fluctuating flux liquid phase. 

Columnar defects too may be naturally present in the 
sample as line dislocations, and can in addition be artifi¬ 
cially introduced by growing them using, for instance, 
MgO nanorods |I7j | or through irradiating the sample 
with high-energy ions such as Sn, Pb, or Iodine. In 
contrast with point-like disorder, the presence of corre¬ 
lated columnar defects results in the emergence of a novel 
low-temperature thermodynamic state distinct from the 
vortex glass p hase, namely the strongly pinned Bose 
glass [3, IlSl - l^ . Since vortex lines become then localized 
along the entire length of these linear pinning centers, the 
sample’s tilt modulus diverges (transverse Meissner ef¬ 
fect) [iO, [13' Indeed, columnar defects have experimen¬ 
tally proven more efficient at pinning than uncorrelated 
point-like disorder due to their extended nature [^ . 

We emphasize that in the vortex glass phase in the 
presence of point-like disorder, the low-energy flux line 
contours tend to roughen as the vortices minimize their 
free energy in the disordered pinning landscape, while in 
the Bose glass phase in the presence of columnar defects 
that are oriented along the magnetic field direction the 
bound vortex lines are effectively straight. Therefore, 
these two disorder-dominated glassy phases in fact dis¬ 
play very different lateral and transverse vortex line fluc¬ 
tuations. One consequently expects samples with point 
disorder and columnar defects to yield profoundly dis¬ 
tinct non-equilibrium relaxation properties. 

Some materials undergo slow relaxation processes to 
eventually reach thermal equilibrium, which affects their 
measured dynamical properties and renders these effec¬ 
tively time-dependent; this is referred to as physical ‘ag¬ 
ing’ m, |2^. Many glassy systems have been observed 
to undergo physical aging [26| owing to the very slow 
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kinetics between a large number of energetically close 
metastable states typically emerging in frustrated envi¬ 
ronments. In a superconducting 2H-NbSe2 sample, Du 
et al. discovered that the voltage response to an ap¬ 
plied current pulse depended on the pulse duration [27|, 
which constitutes clear evidence of physical aging in dis¬ 
ordered vortex matter. Papadopoulou et al. observed 
aging effects in the measurement of the zero-field cooled 
magnetization in the compound Bi 2 Sr 2 CaCu 208 [^ . 

On the numerical front, Bustingorry, Cugliandolo, and 
Dominguez pioneered Langevin molecular dynamics sim¬ 
ulations of a three-dimensional model of vortex matter to 
identify physical aging fe atures in certain two-time cor¬ 
relation functions |29l . l30| . Langevin molecular dynamics 
simulations were also performed by Cao et al. to study 
the dynamical scaling properties of a three-dimensional 
vortex line model subject to an ac drive near the criti¬ 
cal depinning threshold [Slj. Pleimling and Tauber then 
employed an elastic line model and Monte Carlo simula¬ 
tions to investigate the non-equilibrium relaxation prop¬ 
erties of magnetic vortices in type-II superconductors 
starting from somewhat artificial initial configurations 
where straight flux lines were placed at random posi¬ 
tions in the sample [s^. The resulting complex aging 
features (with identical initial conditions and parameter 
values) were subsequently confirmed in a very different 
microscopic representation of the non-equilibrium vortex 
kinetics through Langevin molecular dynamics [s^. A 
systematic investigation of the rich non-equilibrium re¬ 
laxation dynamics and aging features of complex systems 
such as disordered vortex matter with various competing 
energy and time scales is of fundamental interest. Specif¬ 
ically, it is crucial to assess the relevance of distinct initial 
configurations and their effects on the subsequent slow 
relaxation kinetics. 

In this present work, we investigate in detail the ag¬ 
ing relaxation dynamics of vortex matter in disordered 
type-II superconductors starting from initial conditions 
that more closely resemble experimentally realizable sit¬ 
uations, where we analyze the effects of sudden changes 
in either the system’s ambient temperature or magnetic 
field. This allows us to both make closer contact to ex¬ 
periment, and to explicitly test the sensitivity of the non¬ 
equilibrium aging kinetics to the selection of initial states. 
As in our earlier studies, we can switch off the repulsive 
vortex interactions, the pinning potentials, or both in our 
computer simulations, and thereby identify the physical 
mechanisms for the observed complex features in the vor¬ 
tex relaxation dynamics. Again, crucial differences in the 
fluctuation spectrum for vortex matter in the presence of 
point-like or extended defects become apparent. In ad¬ 
dition, separately tracking the behavior of the originally 
present, already relaxed flux lines and the newly added 
vortices helps our analysis of the simulations with sudden 
magnetic field increases. 

One of the principal goals of this investigation is to 
identify observable correlations that allow a dynami¬ 
cal characterization of material samples. Experimental 


methods that could be used to extract the relevant quan¬ 
tities include the following techniques: Vortices in type-II 
superconductors were imaged by Vasyukov et al. using 
nano-SQUIDs in scanning probe microscopy [s^, with 
recorded magnetic fields down to 50 nT. Furthermore, 
Auslaender et al. used the tip of a magnetic force micro¬ 
scope to drag the end of a single vortex across the sur¬ 
face of a YBCO sample, directly measuring its interaction 
with a local disorder potential, in order to study the dy¬ 
namics of vortices and their (de-)pinning processes from 
different defects [s^. Flux creep parameters in YBCO 
samples were experimentally investigated by Abulafia et 
al. by utilizing an array of microscopic Hall sensors . 
Structural properties of vortices in superconductors can 
be imaged through small-angle neutron scattering (3^ . a 
technique that directly measures the Fourier transform 
of the flux lines’ height-height correlation function, and 
thus accesses their lateral fluctuations reflected in the 
observables introduced in Section III El These methods 
could perhaps be utilized to probe the relaxation pro¬ 
cesses of magnetic vortex lines in a disordered medium 
following quenches in temperature or magnetic field, and 
to distinguish the effects of different types of disorder on 
the dynamics of this complex system. The initial con¬ 
figuration of randomly placed straight vortices that were 
implemented in earlier numerical studies cannot be real¬ 
ized in experiments. Instead, we investigate the effects 
of abrupt changes of external experimental control pa¬ 
rameters, namely the system’s ambient temperature and 
magnetic field. 


This paper is organized as follows: In the next sec¬ 
tion, we describe the elastic line model and the Langevin 
molecular dynamics simulation algorithm that we em¬ 
ploy in our study, and also specify the material parame¬ 
ters that we use in our computational model. Then, we 
specify the single-time quantities and two-time correla¬ 
tion functions that are measured in our simulations to 
analyze the non-equilibrium and aging properties of our 
vortex matter models. The last part of this section is 
devoted to a discussion of finite-size effects and a verifi¬ 
cation that our numerical results are not specific to the 
system size we choose. Section IhT] discusses the effects 
of temperature quenches on the relaxation properties in 
systems of interacting flux lines subject to either point¬ 
like or columnar disorder. In Section lYl we present the 
relaxation dynamics in the aging regime following mag¬ 
netic field quenches. We compare the effects of sudden 
increases or decreases of the vortex density to our ear¬ 
lier studies where the magnetic field remained fixed. We 
systematically disentangle the contributions due to repul¬ 
sive vortex interactions and pinning to disorder. We find 
the non-equilibrium relaxation features of vortex matter 
to strongly depend on the type of disorder present in the 
system. We finally summarize our work in the concluding 
Section |Vl 
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II. MODEL DESCRIPTION AND SIMULATION 
PROTOCOL 


B. Langevin Molecular Dynamics 


A. Elastic Line Description 


In the following numerical study, we consider systems 
composed of N vortex lines in the extreme London limit, 
where the London penetration depth A is much larger 
than the coherence length which allows us to employ a 
fully three-dimensional elastic line description 1 ^ . 
The vortices are described through their trajectories 
ri(z), i.e., their two-dimensional ab plane position vec¬ 
tors as function of the z coordinate along the direction of 
the applied external magnetic field. The system’s Hamil¬ 
tonian can then be written as a functional of the N in¬ 
teracting trajectories fi{z): 


i=l ' 


f dz 

ei 

dr*(z) 

Jo 

2 

dz 


+ UD{ri{z),z) 


1 ^ 

+ {\n{z) - rj{z)\) 




( 1 ) 


This effective free energy consists of three compet¬ 
ing terms: the elastic vortex line tension, an attrac¬ 
tive potential describing the Nd pinning centers, and 
the repulsive interactions between different vortices. The 
elastic line stiffness or local tilt modulus is given by 
ei Ri r“2eoln(Aa6/Ca6), where T”^ = Mab/M^ is the 
anisotropy parameter (effective mass ratio). Xab and ^ab 
respectively denote the London penetration depth and 
the coherence length in the ab crystallographic plane. 
Moreover, the in-plane vortex repulsive interactions are 
given by V{r) = 2eoKo{r/Xab), with Kq denoting the 
zeroth-order modified Bessel function, which essentially 
represents a logarithmic repulsion that is exponentially 
screened at the scale A. We set a cut-off for the interac¬ 
tions in our simulations at 5Xab to avoid artifacts due to 
the periodic boundary conditions. 

We discretize the Hamiltonian © in the z direction, 
and model N]j point-like pinning centers as spatially ran¬ 
domly distributed smooth potential wells of the form 

Nd , 

UD{r,z) = -'^-^pS{z - Za) 

( 2 ) 

where Tq, and Za respectively denote the in-plane and z 
coordinates of pinning center a , and p > 0 represents 
the pinning potential strength. Columnar defects are set 
up by first randomly selecting the in-plane coordinates 
rA, and then placing identical potential wells dD) at all z 
positions z^. Henceforth, all lengths are measured rela¬ 
tive to the pinning potential width fog, and energies are 
measured in units of egbg with eg = (i^g/dTrAab)^, where 
d’g = /ic/2e is the magnetic flux quantum. 


1 —tanh I 5 


’ - - bi 


We analyze the vortex kinetics by means of a Langevin 
molecular dynamics algorithm. As mentioned above, we 
discretize space in the z direction (along the external 
magnetic field) into layers, whose spacing corresponds to 
the unit cell size Cg along the crystal’s c direction [sl, [s^ . 
The effective forces acting on each flux line element are 
then computed as spatial derivatives of the discretized 
version of the Hamiltonian ©■ Correspondingly, we dis¬ 
cretize and numerically solve the overdamped Langevin 
equation 


dn{t,z) 


5H[ri{t,z)\ 

Sri{t,z) 


+ Mt,z ), 


( 3 ) 


with the Bardeen-Stephen viscous drag parameter t] by 
direct temporal integration [i^. Note that all fast de¬ 
grees of freedom (e.p., phonons etc.) are captured by 
thermal stochastic forces, which are modeled as un¬ 
correlated Gaussian white noise with vanishing mean 
(fi(t,z)) = 0. We impose the Einstein relation 


z)fjis, z')^ = 2r]kBT5ij5{t - s)(5(z - z') (4) 

to guarantee that the system asymptotically reaches ther¬ 
mal equilibrium with a canonical probability distribu¬ 
tion oc ^-HlksT^ lines are moving in a three- 

dimensional space with periodic boundary conditions in 
the xy plane and free boundary conditions in the z direc¬ 
tion. We set the system size to (16/-\/3)Aah x 8Aah, since 
this aspect ratio ensures that the system equilibrates to 
form a hexagonal Abrikosov lattice of the present flux 
lines in the absence of disorder and driving current. 


C. Material Parameters 

We select our material parameters guided by the ac¬ 
tual values for the ceramic high-Tc type-H superconduct¬ 
ing compound YBa 2 Cu 307 (YBCO). All simulation dis¬ 
tances are measured in units of the pinning center radius, 
which we set to 5g = SSA. We choose the spacing between 
layers in the c direction equal to this basic microscopic 
scale, Cg = 6g. The in-plane London penetration depth 
is Xab = 346g r: 1200A, and the superconducting coher¬ 
ence length ^ab = 0.36g r 10.5A. This emulates the high 
anisotropy of YBCO, with an effective mass anisotropy 
ratio = 1/5. All simulation energies are measured in 
units of eg5g, where eg is the line energy per unit length: 
eg R 1.92 • 10“®erg/cm. The vortex line tension en¬ 
ergy scale becomes ei/eg R 0.189. The depth of the 
pinning potential wells is fixed at p/eg = 0.05. We per¬ 
form our numerical runs at temperatures T = 10 ... 40 K, 
with T = lOK corresponding to kuT/eobo = 0.002 
in our simulation units. Finally, we set the Bardeen- 
Stephen viscous drag coefficient y = 4>o/‘^T^Pnc'^iab ~ 












4 


10 ^°erg-s/cm^ to one (for the normal-state resistiv¬ 
ity of YBCO near T^, Pn ~ 500/^17 cm, see Table 1 in 
Ref. ill); thus the simulation time steps are fixed by 
the fundamental temporal unit to = b^o/eo ~ 18 ps. In 
the following, all times are measured in simulation time 
steps, i.e., in units of to- 


D. Simulation Protocol 


Our aim is to impose experimentally realizable initial 
conditions to bridge theoretical simulations and experi¬ 
mental work. To this end, we investigate sudden changes 
in the magnetic field or temperature. A magnetic field 
quench is implemented in our simulations by adding flux 
lines at random positions in the system (up-quench) or 
removing randomly selected vortices from the system 
(down-quench). We start with straight flux lines posi¬ 
tioned at random in the sample, and let the system re¬ 
lax up to an initial relaxation time r = 10®. We then 
instantaneously increase the magnetic field by adding 5 
straight flux lines to the 16 originally present and essen¬ 
tially relaxed vortices; or instantaneously decrease the 
field by removing 5 random lines from 21 relaxed initial 
vortices. Subsequently, we let the resulting system relax 
up to a waiting time s, at which a snapshot of the sys¬ 
tem is taken; the waiting time s is measured after the 
quench. At later times t > s, different two-time correla¬ 
tion functions are evaluated. Since r is much larger than 
the different waiting times s and measurement time t, it 
is crucial to take this elapsed time period into account 
when considering the initial lines that have relaxed for a 
duration r. We therefore introduce the times cr = r -|- s 
and r = r + t. 

When the external magnetic field is altered in experi¬ 
mental samples, flux lines enter or leave the sample from 
its boundaries. Since our numerical study is limited to a 
fairly small number of vortices, introducing physical sur¬ 
faces through the removal of our implemented periodic 
boundary conditions would render our system dominated 
by the ensuing boundary effects. We rather consider our 
simulation domain as a representative sample of the bulk 
region in a much larger physical system. Also, the forma¬ 
tion or deletion of vortices from the sample’s boundaries 
occurs on time scales much faster than considered in this 
study of comparatively slow relaxation kinetics. 

Temperature quenches are accomplished by letting the 
initial configuration of A = 16 flux lines relax up to a 
similar initial relaxation time (t = 10®) at which the tem¬ 
perature is suddenly raised or lowered. Throughout our 
analysis of magnetic field and temperature quenches, the 
number of layers is L = 640 and the number of pinning 
sites per layer is fixed at Nd/L = 1710. 


E. Measured Quantities 

We measure various two-time correlation functions, 
some of which provide insight into the local features of 
the fluctuating vortex lines, such as their height-height 
autocorrelations 

C{t,s) = ([ri,^(f) - n{t)\ ■ [n,4s) - ri(s)]) , (5) 

where ri,z{i) are the in-plane coordinates of line i at layer 
z at time t, and fi{t) = {ri,z{t)) is the mean lateral po¬ 
sition of line i at time t. Here, observable quantitities 
are averaged over all line elements as well as many noise 
histories and disorder realizations: The brackets (a) in¬ 
dicate the average over noise and disorder realizations, 
where a = Sili denotes the average of 

the local variable Oi^z over all line elements present in 
the system. C{t,s) thus measures the thermal fluctua¬ 
tions of line elements around their corresponding mean 
lateral line position. We remark that the term ‘height- 
height autocorrelation function’ for the transverse vortex 
fluctuations derives from viewing the flux lines as fluctu¬ 
ating one-dimensional interfaces described by a unique lo¬ 
cal height function that measures the deviation of ri,z{t) 
from the respective line’s mean. 

The global structure of the flux line configuration is 
analyzed using the two-time mean-sguare displacement 

B{t, s) = (yi,zit) - ri,z{s)?) ■ (6) 

This correlation function measures the average square 
distance between the position of a line element at time s 
and the position of that same line element at a later time 
t. Both two-time correlation functions were also analyzed 
in Refs. [lH and [s^, where Monte Carlo and Langevin 
molecular dynamics simulations, respectively, were per¬ 
formed to characterize the complex non-equilibrium re¬ 
laxation of initially straight and randomly placed vortex 
lines in the presence of point-like disorder. Thermal spa¬ 
tial fluctuations can also be quantified by means of the 
vortex line radius of gyration 

rg{t) = , (7) 

i.e., the root mean-square displacement from the lines’ 
mean lateral position at time t. 

Another one-time observable that we measure in our 
simulations is the fraction ip{t) of pinned line elements 
when disorder is present. To this end, we first find the 
number of line elements that reside within a fixed small 
cut-off distance from a certain pinning center at a given 
time t. Then we divide this count by the total num¬ 
ber of line elements present in the sample to yield <p(t). 
Throughout our study, the cut-off distance is set equal to 
the pinning center radius, rc = bo. The total number of 
line elements can by computed by multiplying the num¬ 
ber of layers L = 640 in the z direction with the number 
N of vortex lines present in the system. 
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F. Finite-Size Effects 

In our work, we restrict ourselves to a low magnetic 
field or vortex density regime, and only consider systems 
composed of between = 10 and = 21 vortices. This 
naturally constitutes a tiny system compared to typical 
experimental samples of disordered type-II superconduc¬ 
tors in external magnetic fields. It is therefore crucial to 
analyze the effects of this small vortex density, especially 
when we remove lines as the magnetic field sudden de¬ 
creases. We consequently ran simulations in systems with 
different numbers of flux lines, but keeping the density 
of pinning centers the same in all cases. 

As demonstrated in Fig. [1] the data for the vortex ra¬ 
dius of gyration rg{t), mean-square displacement B(J,, 0), 
and the fraction ip{t) of pinned line elements that cor¬ 
respond to systems of different number of flux lines col- 
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FIG. 1: (Color online) Relaxation of the (a) gyration radius, 
(b) mean-square displacement, and (c) fraction of pinned line 
elements for (dilute) systems with different number of inter¬ 
acting vortex lines in the presence of point-like pinning cen¬ 
ters, averaged over 100 disorder and noise realizations. 


lapse onto a common curve. Therefore, we conclude that 
these differently-sized systems follow similar relaxation 
processes confirming that finite-size effects do not play 
a prominent role in this computational study of dilute 
systems of flux lines in a dissipative medium. It is worth 
noting that the physics in denser vortex systems could of 
course be quite different. 


III. TEMPERATURE QUENCHES 

We first consider sudden changes in the ambient tem¬ 
perature of the system of interacting vortex lines, with 
the goal to analyze the effects of instantaneous increases 
in temperature on the system’s relaxation dynamics. At 
t = 10® in our simulation, we instantaneously increase 
the temperature from T = 0.002 to T = 0.008 in a sys¬ 
tem of A^ = 16 interacting flux lines in the presence of 
either point or columnar pinning centers. We then com¬ 
pare the relaxation of the radius of gyration rg{t) and 
the fraction of pinned line elements ip{t) when the tem¬ 
perature stays fixed with the scenario wherein it is sud¬ 
denly raised. We remark that the error bars displayed 
in some of the figures in this section, however small they 
appear, always represent the standard error of the mean 
of the computed quantities. We observe that the relax¬ 
ation of these observables in the considered systems can 
invariably be mapped to an exponential function of the 
form a 6“*/’’’ + c, where t determines a characteristic re¬ 
laxation time. The measured values for the relaxation 
times T in different scenarios are summarized in Table ID 
We note that the exponential relaxation of observables in 
our vortex model system following sudden temperatures 
quenches implies that physical aging features would not 
be discernible on time scales larger than the relaxation 
time T. 


Defect 

type 

T = 0.0021 T = 0.0081 T = 0.002 ^ 0.008 

rg 

‘P 

(a) Point-like 

5.9 • 10^ 

5.3 ■ 10^ 

3.5 ■ 10'’^ 

3.1 ■ 10^ 

(b) Columnar 

3.5 ■ 10^ 

4.9 ■ 10'‘ 

n = 1.6 • 10® 

T2 = 5.7 • 10^ 

3.3 • 10^ 


TABLE I: Measured characteristic relaxation times for the 
gyration radius Vg and the fraction of pinned line elements 
in systems of interacting vortex lines with (a) point-like and 
(b) columnar defects (data averaged over 1000 realizations). 

In systems with point-like disorder, the measured re¬ 
laxation time of the radius of gyration at high temper¬ 
ature is shorter than that at low temperature, see Ta¬ 
ble H This confirms that thermal fluctuations facilitate 
the roughening of vortices in the presence of point-like 
disorder. Following a sudden increase in temperature, 
the radius of gyration first decreases a little as a few 
flux line segments become thermally depinned, and sub¬ 
sequently increases monotonically, as shown in Fig. I2al 
We checked that this increase in rg{t) is well described 
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FIG. 2: (Color online) Relaxation of the gyration radius in 
systems of interacting vortex lines with (a) point disorder or 
(b) columnar defects, when the temperature is either held 
fixed or instantaneously raised (data averaged over 1000 real¬ 
izations). The inset in (a) shows the gyration radius following 
the sudden temperature increase in the presence of point dis¬ 
order, and a fit to an exponential function. 


by an exponential function of the form a -I- c, with 
relaxation time r = 3.5 • 10'^ (the other parameters are 
a ~ —1.73 and c ~ —0.09) as observed in the inset in 
Fig. mi Thus, as expected, point defects enhance ther¬ 
mal wandering of flux lines in this disordered landscape. 

In contrast, systems with columnar defects actually re¬ 
quire a longer relaxation time for the radius of gyration 
at higher temperatures than that at lower temperatures, 
see Table U asserting that thermal fluctuations resist the 
straightening of vortex lines as they become localized at 
columnar pinning centers. When we instantaneously in¬ 
crease the temperature in this system, the gyration ra¬ 
dius after the quench shows a marked fast increase, even 
beyond the relaxation curve for a system at constant 
higher temperature, followed by a much slower decrease. 


see Fig. [2bl The initial increase in Vg can likely be at¬ 
tributed to the creation of double-kinks in the system 
that are not similarly generated in relaxation processes 
at fixed temperature, while the final decrease is due to 
the decay of these metastable configurations. A double¬ 
kink refers to a flux line configuration where the vortex is 
simultaneously pinned to two adjacent columnar defects. 
This is a long-lived state which will eventually decay to 
either a free vortex, or, more probably, to a straight flux 
line bound to a single defect column. We find that both 
the initial increase and the final decrease in the gyration 
radius again fit an exponential of the form a + c, 
but with different relaxation times (and signs of the co¬ 
efficient a): In Fig. I2b[ the fast relaxation time from the 
initial non-equilibrium steady state at T = 0.002 to the 
metastable state with vortex double-kinks is Ti = 1.6T0^, 
while the much longer relaxation time from this interme¬ 
diate state to the final non-equilibrium steady state at 
T = 0.008 is t 2 = 5.7 • 10^. Note that the data for the 
gyration radius for both the temperature up-quench sce¬ 
nario as well as at a fixed high temperature in Figs. 
and I2bl will eventually approach the same curve in the 
long-time limit, which has not yet been reached in our 
simulation time window. 

The differences between the systems’ responses to tem¬ 
perature up-quenches in the presence of point and colum¬ 
nar pinning centers, respectively, can be attributed to 
the very distinct (de-)pinning behavior that we shall 
again discuss later in the context of sudden magnetic 
field changes. When the temperature is suddenly raised 
in systems with uncorrelated point disorder, some line 
elements become thermally depinned. The vortices are 
initially straightened owing to their line tension. Sub¬ 
sequently an increasing number of line elements finds 
favorable pinning sites, causing a monotonic growth in 
their transverse spatial fluctuations measured by Vg. In 
samples with columnar defects, following a temperature 
up-quench some line elements migrate away from certain 
columnar pinning centers, which renders them subject 
to thermal wandering. Yet they soon become pinned to 
neighboring defect lines causing the subsequent decrease 
in the lateral spatial vortex line fluctuations, akin to sys¬ 
tems relaxing at fixed temperature. 

This scenario is confirmed by the relaxation data for 
the fraction of pinned line elements ip{t) in Fig. [3] As 
observed in Fig. I3al in the presence of point disorder 
only a minority (~ 25%) of line elements are actually 
pinned right before the temperature quench. A sudden 
increase in temperature releases a significant fraction of 
these pinned line elements causing a decrease in ip as 
well as the small dip in Vg. Later, the vortices become 
attached to and hence stretched between point-like pins 
causing a subsequent increase in ip, accompanied with a 
monotonic increase in the lines’ spatial fluctuations. 

In the presence of columnar defects, on the other hand. 
Fig. |3b] shows that almost all line elements are bound 
to pinning centers just before the temperature quench 
is applied. The sudden temperature increase causes a 
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FIG. 3: (Color online) Relaxation of the fraction of pinned 
line elements in systems of interacting vortex lines with (a) 
point disorder or (b) columnar defects, when the temperature 
is either held fixed or instantaneously raised (data averaged 
over 1000 realizations). 



FIG. 4: (Color online) Relaxation of the gyration radius in 
a system of interacting vortex lines with columnar pinning 
centers, when the temperature is held fixed or instantaneously 
lowered (data averaged over 1000 realizations). 

perature, where we reduced the system’s ambient tem¬ 
perature at t = 10® from T = 0.008 to T = 0.002 in our 
simulation runs. The subsequent vortex relaxation kinet¬ 
ics displays some noticeable differences as compared with 
the previous up-quenches: Following the temperature de¬ 
crease, the gyration radius decreases and the fraction of 
pinned line elements increases in samples with both types 
of disorder, and non-monotonic behavior is observed in 
neither of these quantities, as depicted for Vg in Fig. ID 
When the temperature suddenly drops, more line ele¬ 
ments simply become pinned to either point-like or ex¬ 
tended defects, which causes a concomitant decrease in 
the transverse spatial fluctuations of the vortex lines. 


IV. MAGNETIC FIELD QUENCHES 


tiny number of these line elements to depin from their 
localizing columns. However, they are soon to be pinned 
again, but partly in metastable double-kink configura¬ 
tions. This explains the non-monotonic pattern we ob¬ 
serve in the spatial fluctuations with a rapid increase 
of Tg followed by a slow decrease. Finally, the fraction 
of pinned line elements ip{t) following the temperature 
quench in the presence of both types of disorder once 
again fits an exponential of the form a + c, yet with 
different relaxation times: r = 3.1-10’^ for point disorder, 
whereas r = 3.3 • 10"* for columnar defects, see Table ID 
This confirms that enhanced thermal fluctuations at el¬ 
evated temperatures resist the straightening of flux lines 
as they attempt to bind to columnar defects; in stark 
contrast, thermal effects facilitate line roughening when 
vortices become pinned to uncorrelated point-like disor¬ 
der. 

We also analyzed the effects of a sudden drop in tem- 


A. Non-interacting Vortex Lines without Disorder 


We now proceed to study the effects of sudden changes 
in the external magnetic field, i.e., vortex density, while 
keeping the ambient temperature fixed at T — 0.002 in 
all simulation runs for the various systems described in 
this section. The dynamics of non-interacting directed 
lines in the absence of disorder can be directly mapped to 
the one-dimensional Edwards-Wilkinson interface growth 
model (which is in turn equivalent to a free noisy diffusion 
equation). The two-time height-height autocorrelation 
function ([5]) in the correlated regime is found to be [i^ 


C(t,s) = Cos^/^ 




( 8 ) 


Indeed, if a simple aging scenario applies, one expects the 
general scaling form for the height-height autocorrelation 





































FIG. 5: (Color online) Relaxation of (a) the mean-square dis¬ 
placement and (b) the height-height autocorrelation function 
for non-interacting vortices in the absence of disorder and for 
fixed vortex density (data averaged over 800 realizations). 

function 

C{t,s) = s'^fc{t/s) (9) 

with a scaling function fc- In the free diffusive Edwards- 
Wilkinson regime, the scaling exponent is & = 1/2 [32], 
Id^ . In our Langevin molecular dynamics study, we first 
analyze this case of free non-interacting vortex lines and 
start with the fixed magnetic field scenario to validate 
our numerical code. Since the present flux lines relax for 
a long initial time r = 10®, which is much larger than 
the waiting times s and measurement time t, we need to 
take this elapsed time r into account when analyzing the 
behavior of these lines, and therefore new times a = r + s 
and r = r-|-t are introduced. The two-time mean-square 
displacement (|6]) follows a similar scaling form to ([9]). As 
shown in Fig. [51 our data for both i?(r, a) and (^(r, a) 
indeed satisfy dynamical scaling for all explored waiting 
times with the predicted scaling exponent b = 1/2, after 


properly accounting for the long initial relaxation time 
r, see Sec. iUDi 

When we follow our magnetic field quench scenario, 
and suddenly reduce the vortex density after letting the 
system relax for an extended time period, we obtain 
the very same relaxation results for both B(r, a) and 
(^(r,^), which is to be expected since the lines are not 
interacting. 

Similarly, when we instantaneously increase the mag¬ 
netic field, B(T,a) and C{r,a) measured for those vor¬ 
tices that were present initially are not significantly dif¬ 
ferent from those quantities obtained at fixed flux den¬ 
sity, as the relaxation processes for these non-interacting 
vortex lines are statistically the same in both situations. 
However, the subpopulation of newly added lines displays 
a similar, but less perfect collapse for C{t,s) than that 
displayed in Fig. |5b| since these lines were not present 
for the previous long relaxation time r, and thus are not 
yet as relaxed as the initial lines. On the other hand, 
B{t,s) obtained just for the added lines exhibits perfect 
dynamical scaling with the scaling exponent b = 1/2. 

In summary, we confirm the expected Edwards-Wilkin- 
son aging scaling for the non-equilibrium relaxation dy¬ 
namics of non-interacting flux lines in the absence of dis¬ 
order, regardless of sudden changes in their density. The 
following subsections highlight the effects of disorder on 
the relaxation dynamics of non-interacting, thus statisti¬ 
cally independent vortex lines. 

B. Non-interacting Flux Lines with Point Disorder 

Keeping the repulsive vortex-vortex interactions 
switched off, we now introduce point defects into the sys¬ 
tem to study the effects of this type of spatially uncorre¬ 
lated disorder on the relaxation dynamics. 



FIG. 6: (Color online) Relaxation of the normalized height- 
height autocorrelation function in a system of non-interacting 
flux lines subject to point-like disorder, when the number of 
flux lines stays fixed (data averaged over 800 realizations). 
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First, we consider the situation of fixed magnetic field 
(vortex density), and study the same dynamical quanti- 
ties as in the previous section, see also Refs. [^[^. The 
mean-square displacement i?(r, tr) does not display dy¬ 
namical scaling with & = 1/2 as in the disorder-free case. 
Instead, dynamical scaling ensues for i?(r, tr) with the 
aging exponent b r; 0.725, c.f. Fig. 6(d) in Ref. [s^. In 
contrast, dynamical scaling is not observed for the height- 
height autocorrelation function ^(F, cr). As shown in 
Fig. [HI the data for all waiting times s display identical 
time evolution, since the system has in fact fully equi¬ 
librated owing to the long relaxation period r = 10®, 
and time translation invariance is restored. The height 
autocorrelation function shows a similar increase for all 
waiting times indicating that flux line elements keep ex¬ 
ploring their surroundings in search of favorable pinning 
sites, which causes the lines to continuously roughen due 
to the uncorrelated nature of the point-like disorder in 
the system. 

The results for a sudden magnetic field decrease (vor¬ 
tex density down-quench) turn out to not significantly 
differ from the situation with constant field. This is read¬ 
ily attributed to the fact that the lines are not interact¬ 
ing, and thus the remaining vortices after the quench are 
not at all affected by the removal of some lines from the 
sample. 

More distinct features are observed upon an instan¬ 
taneous increase of the magnetic field: First of all, the 
mean-square displacement i?(F, a) of the initial lines does 
not display dynamical scaling anymore, but the data 
(not shown) reveal that this vortex subpopulation has 
rather recovered time translation invariance, since they 
all collapse onto a single curve when plotted against 
F —tr = t — s, while the added line subpopulation displays 
neither of these straightforward features. The height- 
height autocorrelation function (^(F, a) of the initial lines 
shows time translation invariance, akin to the data in 
Fig.ini On the other hand, C{t, s) for the added lines, de¬ 
picted in Fig. [Tal displays an initial decrease followed by 
an increase, with the minimum shifted to later t—s values 
for longer waiting times s. The decrease happens when 
the added lines start noticing the presence of the pin¬ 
ning centers and become pinned to nearby defects, which 
decreases their transverse fluctuations. The subsequent 
increase is caused by the flux lines exploring more of the 
sample and thus becoming rougher as they optimize their 
pinning configurations. 

Hence, the behavior of the added lines is highly influ¬ 
enced by disorder as demonstrated by their height auto¬ 
correlations. A one-time observable that also highlights 
the effects of the sudden increase in the magnetic field is 
the fraction of pinned vortex line elements tp. In the pres¬ 
ence of point disorder, this fraction as shown in Fig. Eb] 
displays a two-step increase, which can be attributed to 
the pinning dynamics of the newly added magnetic flux 
lines, as discussed in the context of Fig. [7al 




FIG. 7: (Color online) Relaxation of (a) the normalized 
height-height autocorrelation function of the newly added flux 
lines, and (b) the fraction of pinned line elements in a system 
of non-interacting vortices subject to point disorder, follow¬ 
ing a sudden increase in the magnetic field or vortex density 
(data averaged over 800 realizations). 


C. Non-interacting Vortices with Columnar Pins 

We proceed to change the type of pinning centers 
present in the system to columnar defects, which are lin¬ 
early extended pinning sites that are simulated as po¬ 
tential wells aligned along the z direction, and compare 
their effects on the relaxation kinetics of non-interacting 
vortices to that of uncorrelated point disorder. 

We again first hold the magnetic field constant and 
consider the different two-time correlation functions 
studied above. The mean-square displacement B(T,a) is 
found to display dynamical scaling with the aging expo¬ 
nent b = 1/2, similar to Fig. [5^ in the absence of disorder. 
It is worth recalling that when flux lines are completely 
pinned to columnar defects, their mean-square displace¬ 
ment becomes constant, of the order of bg. We therefore 
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FIG. 8: (Color online) Relaxation of various observables in a system of non-interacting flux lines with columnar defects, (a,c) 
when the vortex density remains fixed; (b,d) after a sudden increase in the magnetic field: (a,b) normalized height-height 
autocorrelation function, (c,d) fraction of pinned line elements (data averaged over 800 realizations). 


conclude that the dynamics of B{T,a) in the accessible 
time window is clearly dominated by the still unattached 
lines in the system. The height-height autocorrelation 
function (^(r, cr) does not display dynamical scaling with 
any fixed scaling exponent. However, when considering 
the normalized C(r, cr) with respect to the initial value 
C{a,a) as depicted in Fig. Eil we observe that the data 
for all different waiting times s perfectly collapse onto 
a single curve, which indicates that time translation in¬ 
variance is restored and the system is well equilibrated. 
Next, we measure the fraction of pinned line elements 
in the system with columnar defects, see Fig. [53 The 
data for different waiting times collapse onto a single 
graph, where the fraction of pinned line elements stays 
constant for an extended time period and only later, at 
ln(r — cr) Ri 7, starts increasing. This indicates that 
F Ri 1096 + cr represents a characteristic time scale when 
pinning to columnar defects becomes prominent in the 
system. 


When vortices are removed from this system of non-in¬ 
teracting flux lines with columnar defects, the measured 
mean-square displacement, height autocorrelations, and 
fraction of pinned line elements all yield identical results 
to the situation when there is no sudden magnetic field 
change, as is to be expected in the absence of correlations. 

Upon introducing additional flux lines into the system, 
the mean-square displacement of both the initial and the 
added lines subpopulations shows dynamical scaling for 
all waiting times with the same exponent 5 = 1/2 as in 
the fixed density and field down-quench scenarios. This 
confirms that sudden changes in the vortex density do not 
alter the global dynamics of non-interacting flux lines in 
the presence of columnar defects. The autocorrelation 
function C(r, cr) for the initially present lines turns out 
very similar to that at fixed density, see Fig. [53 whereas 
the data corresponding to the newly added lines cannot 
be collapsed according to a simple aging scaling form 
when plotted against the time ratio t/s. However, the 
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graphs C(t, s) for the added lines for various s do not 
just fall onto a single curve when plotted against t — s, 
as demonstrated in Fig. |8bl This indicates that this sub¬ 
population of lines is not yet equilibrated for the range 
of waiting times s employed here. As the additional lines 
are inserted into the sample, depending on the elapsed 
waiting time s, they require different time periods t to ex¬ 
plore the system and become attached to columnar pins. 
This is corroborated in the time evolution of the fraction 
of pinned line elements cp for different waiting times in 
Fig. |8dl which also does not yield data collapse. 

The flux lines present in the sample, especially vortices 
that were suddenly added to the system, start exploring 
the surrounding pinning centers sooner for shorter wait¬ 
ing times, whence one observes an earlier increase in the 
fraction of pinned line elements for smaller s. For the long 
waiting times in Fig. I8dl such as s = 2^^, the fraction of 
pinned line elements is fixed at a plateau value for an ex¬ 
tended duration because the vortices already had enough 
time to relax towards favorable configurations after the 
field up-quench. The increase in ip{t) happens at a much 
later time indicating that the line elements might be mi¬ 
grating between neighboring columnar pinning centers 
to become pinned to whole defect lines, whereupon their 
lateral line fluctuations diminish. The overall trend in 
the presence of columnar defects is that the dynamics is 
dominated by the few lines that remain unattached to 
any pins, and the tendency of each flux line to become 
fully localized to a single columnar pinning center, which 
in turn straightens the vortex lines. 

The preceding discussions highlight the influence of dif¬ 
ferent disorder on non-interacting vortex matter and the 
relaxation dynamics following magnetic field quenches. 
We expect that quite distinct relaxation kinetics will be 
observed when mutual vortex interactions are added to 
our consideration and collective behavior ensues, as will 
be explored in the following subsections. 


D. Interacting Vortex Lines without Disorder 

We start with the system of interacting flux lines with¬ 
out disorder. At a fixed magnetic field, we observe at 
least approximate dynamical scaling of B{T,a) as dis¬ 
played in Fig. 1211 yet with a different exponent b = 1.57 
from before. On the other hand, the height-height auto¬ 
correlation function in Fig. I9bl shows a distinct two-step 
relaxation behavior, where a /3 relaxation with almost no 
changes in C(F, tr) is followed by a very slow decay. This 
is referred to as a-P relaxation scenario, a characteristic 
feature of structural glasses [11,143. Furthermore, a local 
maximum in C(F, cr) is observed for all waiting times s, 
which we attribute to the decay of metastable kink con¬ 
figurations. Such single kinks are step-like flux line struc¬ 
tures that are stabilized by the screened logarithmic re¬ 
pulsive forces; hence they constitute interaction-induced 
features that are long-lived. The mean-square displace¬ 
ment B{T,a) loses the dynamical aging scaling property 




FIG. 9: (Color online) Relaxation of (a) the mean-square dis¬ 
placement and (b) the normalized height-height autocorrela¬ 
tion function in a system of interacting flux lines in the ab¬ 
sence of disorder, when the number of flux lines stays fixed 
(data averaged over 800 realizations). 


when lines are removed from the system, whereas (^(F, a) 
displays similar behavior to that at fixed magnetic field, 
but with the local maxima absent for most waiting times. 
This indicates a destabilization of single-kink structures 
due to the decrease in the number of lines in the sam¬ 
ple, which relieves the effective mutual cages generated 
by the repulsive interactions. 

On the other hand, when the magnetic field is sud¬ 
denly increased and vortices are added to the system, we 
observe the same time evolution in the mean-square dis¬ 
placement of the added lines, the initial lines, and all the 
lines combined. The data corresponding to all waiting 
times collapse onto a single curve with a long-time al¬ 
gebraic growth with the aging scaling exponent b — 0.1, 
see Fig. llOal This dynamical scaling is manifest only 
when the initial relaxation time r is disregarded in our 
scaling formula for both the initial lines and the collec¬ 
tion of all lines. The independence of these subpopu- 
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FIG. 10: (Color online) Relaxation of (a) the mean-square 
displacement of all the flux lines combined, the normal¬ 
ized height-height autocorrelation function of (b) the initially 
present vortices, and (c) the newly added lines in a system of 
interacting flux lines in the absence of disorder, following a 
sudden increase in the magnetic field (data averaged over 900 
realizations). 


lations from the initial relaxation time implies that the 
respective lines do not retain memory of their early dy¬ 
namics and only experience changes upon the addition 
of new vortices, reflecting the mutual interactions and 
consequent dynamical correlations. 

We observe a more diverse temporal evolution in the 
height-height autocorrelation functions. The initial lines 
in Fig. llObI show a-P relaxation for all waiting times s 
with the presence of local maxima similar to those in 
Fig. I9bl The data display a rather complicated depen¬ 
dence on the sequence of waiting times, which we cannot 
explain in detail but attribute to the rich and complex 
nature of the combined vortex system when these initially 
present and largely relaxed flux lines interact with newly 
added straight ones. Unlike the initial line subpopula¬ 
tion, the added vortex kinetics presented in Fig. dis¬ 
plays a plateau followed by an increase, which is in turn 
followed by a final decrease for all waiting times. The first 
increase indicates the time scale when the added lines 
start noticing the repulsive effect of the already existing 
vortices, whence lateral fluctuations are enhanced. At 
later times following this initial increase, C(t, s) decreases 
again indicating that the lines have reached a steady state 
with well-established mutual interaction cages, and hence 
their fluctuations become reduced. 

In the following subsections, we investigate and ana¬ 
lyze the effects of disorder on samples of interacting vor¬ 
tex lines. We first address uncorrelated point-like defects 
and then proceed to extended columnar pins. 


E. Interacting Flux Lines with Point Disorder 

Considering first the case when the magnetic field does 
not suddenly change, our sample of interacting lines in 
the presence of point disorder shows dynamical scaling for 
B{T,a) with the aging exponent b = 1.54, similar to the 
scaling behavior displayed in Fig. in the disorder-free 
interacting vortex system. This similarity suggests that 
the strong repulsive forces between the flux lines domi¬ 
nate over the effects of point-like pinning centers, which 
do not play a significant role in modifying the global re¬ 
laxation dynamics of interacting vortices. Even though 
the height-height autocorrelation function does not show 
dynamical scaling, the data for the different waiting times 
collapse onto a single curve indicating that the system 
has effectively equilibrated, see Fig. Illal The fraction 
of the pinned line elements (p{t) depicted in Fig. Illci too 
shows evidence for equilibration for all waiting times con¬ 
sidered in our simulations. The data display a plateau 
followed by an increase in the pinned line element fraction 
at ln(F —tr) Ri 7, which coincides with the time of marked 
decrease in C{T,a), confirming that F r: 1096 -I- a rep¬ 
resents a characteristic temporal scale when pinning to 
point-like defects causes significant flux line roughening. 

When flux lines are suddenly removed from the system, 
B(T, a) of the remaining lines does not display dynamical 
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FIG. 11: (Color online) Relsixation of various observables in a system of interacting flux lines with point-like disorder when (a,c) 
the vortex density remains fixed, and (b,d) after a sudden magnetic field decrease: (a,b) normalized height-height autocorrelation 
function, (c,d) fraction of pinned line elements (data averaged over 800 realizations. 


scaling. Moreover, C(r, a) in Fig. lllbi still shows the a-/? 
relaxation scenario for all waiting times without the sys¬ 
tem ever becoming fully equilibrated as in Fig. Illal The 
fraction of pinned line elements undergoes a more diverse 
temporal evolution for the different waiting times, as can 
be observed in Fig. Illdl For short waiting times s < 2®, 
a minimum in ip occurs at ln(r — cr) ~ 6.2. The initial 
decrease indicates that the sudden release of caging due 
to removed lines leads to depinning of some remaining 
vortices from point defects where the line elements were 
only held in place because of their interactions with the 
removed lines. This enables the remaining flux lines to 
explore a farther range of their surroundings, ultimately 
becoming bound to new and more point pins to eventu¬ 
ally reach a steady-state fraction p that is a little larger 
than the initial one. For intermediate-to-larger waiting 
times s > 2®, this minimum in the pinned line element 
fraction is not visible and the data looks similar to the 
situation with fixed vortex density, see Fig. Illcl the only 


difference being that the quenched system has clearly not 
equilibrated yet. 

Let us now explore the effects of a sudden increase 
in the magnetic held, where new straight hux lines are 
added to the a previously relaxed system of interacting 
vortices in the presence of point disorder, and analyze 
the resulting dynamics. Dynamical scaling of the mean- 
square displacement is lost for the newly added lines, the 
initially present vortices, and all the hux lines combined. 
However, the mean-square displacement as function of 
time looks similar among all these different subpopula¬ 
tions; an example is shown in Fig. I12al Note that our 
choice of a putative aging scaling exponent b ~ 0.8 in 
Fig. I12al is merely to demonstrate the global dynamical 
evolution of B{t,s) for the different vortex subpopula¬ 
tions. 

On the other hand, the height-height autocorrelation 
function shows similar results for the initial lines and 
the combination of all vortices, as shown in Fig. I12b[ 
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FIG. 12: (Color online) Relaxation of (a) the mean-square 
displacement, the normalized height-height autocorrelation 
function (b) of all the flux lines combined, and (c) of the 
subpopulation of newly added lines in a system of interacting 
vortices subject to point disorder, following a sudden increase 
in the magnetic field (data averaged over 800 realizations). 


where a-P relaxation is evident as a characteristic fea¬ 
ture in these flux line subpopulations that have already 
relaxed for a long time period, but are still not equili¬ 
brated. This lack of equilibration becomes manifest in 
the height-height autocorrelation function for just the 
subpopulation of added lines in Fig. ll2cl For these newly 
introduced vortices, only very long waiting times s > 2^^ 
display a-P relaxation kinetics comparable to the situa¬ 
tion at fixed magnetic field. However, data corresponding 
to s < 2^^ show an initial increase in C{t, s) followed by 
a slow decrease. When new flux lines are inserted into 
the sample, the originally present vortices repel them and 
thus facilitate their pinning to the uncorrelated point dis¬ 
order, which enhances transverse line fluctuations. At 
ln(t — s) Ri 7, the added lines’ height autocorrelation 
function decreases because the lines have found a (near) 
optimal pinning configuration. The fraction of pinned 
line elements behaves analogously to that in the case 
of a sudden vortex density decrease. Fig. Illdl with the 
mere difference that the minimum is less pronounced for 
a sudden magnetic field up-quench. Yet in this case, the 
sudden presence of additional caging potentials (rather 
than their release) and subsequent vortex rearrangement 
leads to the observed behavior of the pinned line element 
fraction. 


F. Interacting Vortices with Columnar Defects 

Finally, we analyze a sample of interacting flux lines 
in the presence of columnar defects, and compare the 
relaxation kinetics in this system with the preceding in¬ 
vestigation of uncorrelated point disorder. 

We first study the case when the magnetic field stays 
constant. We find that B{t, s) shows dynamical scaling 
at long times with the aging scaling exponent b — 0.8, see 
Fig. 1131 In contrast, no dynamical scaling is obtained for 
C(r, O’) which displays a time evolution similar to that 
in Fig. Illal for point defects. The fraction of pinned line 
elements in the presence of columnar defects turns out to 
not significantly differ from that in the presence of point 
disorder at fixed vortex density. 

When flux lines are removed from the sample, dynami¬ 
cal scaling is not observed for B{T, a), even at long times. 
Furthermore, time translation invariance in the normal¬ 
ized height-height autocorrelation function indicates an 
essentially equilibrated system that undergoes a-P relax¬ 
ation. This behavior is accompanied by an increasing 
fraction of pinned line elements confirming that the re¬ 
maining flux lines undergo similar relaxation dynamics as 
in the system at fixed magnetic field, namely a continuous 
straightening and hence reduced transverse fluctuations 
as the vortices find nearby columnar pinning centers. 

Quite different features are observed when the magne¬ 
tic field is increased instantaneously: The temporal 
growth of the mean-square displacement B(t,s) of the 
newly introduced lines looks similar to samples with 
point-like disorder as shown in Fig. I12al The similar- 
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FIG. 13: (Color online) Relaxation of the mean-square dis¬ 
placement in a system of interacting flux lines with columnar 
defects, when the vortex density remains hxed (data averaged 
over 1500 realizations). 


ity in the mean-square displacement in the presence of 
both point and columnar pinning sites underscores our 
earlier observation that introducing disorder into the in¬ 
teracting vortex system has comparatively little effect on 
the global dynamics which is highly dominated by the 
strong repulsive forces between the flux lines. Analogous 
to the case when point disorder is present, all the dif¬ 
ferent vortex subpopulations show similar results for the 
mean-square displacement. 

On the other hand, the added vortices and the initially 
present flux lines are characterized by different height au¬ 
tocorrelations. The initial vortices as well as the combi¬ 
nation of all lines display a-jS relaxation for all waiting 
times s. The appearance of local maxima for some long 
waiting times indicate the decay of long-lived vortex dou¬ 
ble kinks in the system. However, the subpopulation of 
added lines displays different features in their height au¬ 
tocorrelations C{t,s) for the different waiting times, see 
Fig. I14al For short waiting times, e.g. s = 2®, an initial 
decrease is followed by a pronounced maximum in (7(1, s). 
The initial decrease is set by the time scale when the 
added flux lines begin to notice the pinning centers, caus¬ 
ing their local transverse fluctuations to decrease. The 
following increase indicates the characteristic time range 
when the additional vortices start feeling the repulsive 
interactions with the initially present flux lines and with 
one another; this produces enhanced fluctuations about 
their mean lateral positions. This increase continues un¬ 
til pinning effects start to again dominate the relaxation 
dynamics, rendering the lines straighter and thus reduc¬ 
ing vortex line fluctuations. For longer waiting times, 
a-P relaxation becomes more prominent, with small lo¬ 
cal maxima present for waiting times s = 2® and s = 2®. 
Moreover, the fraction of pinned line elements in Fig. lldbl 
shows an increase followed by asymptotic saturation for 




FIG. 14: (Color online) Relaxation of (a) the normalized 
height-height autocorrelation function of the newly added flux 
lines, and (b) the fraction of pinned line elements in a system 
of interacting vortices with columnar defects, following a sud¬ 
den increase in the magnetic field (data averaged over 800 
realizations). 


all waiting times tested in our simulations. 


V. CONCLUSION AND OUTLOOK 

In this paper, we have explored the effects of incor¬ 
porating experimentally motivated initial conditions into 
the study of the non-equilibrium relaxation dynamics of 
magnetic flux lines in type-II superconductors. Random 
initial conditions were considered in our earlier work in 
Refs. [H, [13 . In the present study, we focused on the 
effects of sudden changes in temperature and magnetic 
field on the ensuing vortex relaxation properties. We in¬ 
vestigated the differences between the relaxation kinetics 
of flux lines in the presence of uncorrelated point-like dis¬ 
order and columnar defects. To this end, we modeled the 
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magnetic vortices as directed elastic lines and employed a 
Langevin molecular dynamics algorithm to simulate the 
vortex kinetics. 

We studied the relaxation of a system of interacting 
flux lines in the presence of disorder at fixed temper¬ 
atures to confirm that higher temperatures induce en¬ 
hanced thermal fluctuations that counteract the straight¬ 
ening of the vortex lines caused by their pinning to colum¬ 
nar defects, while in contrast these same fluctuations help 
lines to roughen when binding to uncorrelated point-like 
disorder. Furthermore, we investigated the relaxation 
processes following sudden changes in the ambient tem¬ 
perature. Systems subject to the two types of disorder 
show very distinct relaxation dynamics when the tem¬ 
perature is instantaneously increased due to the different 
(de-)pinning behavior from these defects. Some flux line 
elements become depinned from the uncorrelated point 
defects following a sudden increase in the sample’s tem¬ 
perature confirming that this type of disorder enhances 
the thermal wandering of flux lines. However, samples 
with columnar defects display a non-monotonic behavior 
in spatial fluctuations due to the creation and subsequent 
decay of vortex double-kinks. These pronounced differ¬ 
ences between systems with point-like and columnar de¬ 
fects are not present when the temperature is suddenly 
reduced, since this decrease simultaneously binds more 
line elements to the pinning centers, regardless of the 
type of disorder present in the sample. 

We carefully investigated the effects of sudden changes 
in the magnetic field, i.e., vortex density, on the relax¬ 
ation processes of flux lines, and performed a systematic 
analysis to disentangle the contributions of the repul¬ 
sive vortex-vortex interactions and pinning to the dis¬ 
tinct point-like and columnar defects. We compared the 
cases when the magnetic field instantaneously increases 
or decreases to the situation when it stays fixed. We start 
with perfectly straight and randomly placed vortex lines 
and let the system relax for an extended relaxation time 
period, after which we maintain a fixed number of lines 
representing a constant magnetic field. Alternatively, we 
remove randomly selected lines from the initially relaxed 
vortices in the sample, or add perfectly straight new flux 
lines at random positions into the system of relaxed lines, 
to analyze the effects of a sudden decrease or increase in 
the magnetic held, respectively. 

We validated our numerical studies by showing that 
our results for free, non-interacting vortex lines agree 
with the analytical predictions from the Edwards- 
Wilkinson interface growth model. We then introduced 
point-like and columnar pinning centers to a system of 
non-interacting hux lines. We found that the relaxation 
features in samples with columnar disorder are quite sim¬ 
ilar to the case when no disorder is present, especially in 
the behavior of global observables, as e.g., visible in the 
dynamical scaling of the vortex mean-square displace¬ 
ment. This similarity indicates that the dynamics of vor¬ 
tex lines in the presence of columnar disorder is strongly 
dominated by the still unattached lines in the system. In 


the three studied cases without vortex interactions, sud¬ 
den removal of flux lines results in the remaining lines 
showing similar relaxation behavior to the situation with 
a fixed vortex density, as expected since caging effects due 
to the repulsive vortex-vortex interactions are absent. On 
the other hand, when additional lines are introduced to 
these systems, the initially existing lines and the added 
lines show quite distinct features in their relaxation dy¬ 
namics. 

Thereafter, we added the repulsive screened logarith¬ 
mic interactions into our consideration of systems with¬ 
out disorder, with uncorrelated point-like, or with ex¬ 
tended columnar defects. Sudden decreases in the mag¬ 
netic field result in the remaining lines displaying relax¬ 
ation properties that differ from the situation when the 
magnetic field stays fixed, which confirms that caging 
effects due to vortex-vortex interactions constitute a 
prominent factor that highly affects the results of sud¬ 
den changes in the magnetic field. Collective interaction 
effects are markedly enhanced when accompanied with 
disorder following sudden increases in the magnetic field, 
where the relaxation properties of the added lines clearly 
depend on the type of pinning centers. Point-like disor¬ 
der enhances transverse spatial vortex line fluctuations 
by binding each to many pinning centers at once in the 
vortex glass phase, while columnar defects pin entire flux 
lines and straighten them in the distinct Bose glass phase. 
It is worth noting that when comparing the novel results 
in this paper to our earlier study in Ref. [s^ with ran¬ 
dom initial conditions, we discovered improved dynami¬ 
cal scaling in the vortex mean-square displacement in all 
situations considered here. Consequently, the currently 
utilized more realistic initial conditions help the system 
to equilibrate faster and thus acquire universal features. 

The goal of this study is to methodically work towards 
observables and protocols that can be realized in exper¬ 
iments, while simultaneously providing detailed theoret¬ 
ical insight into the non-equilibrium relaxation kinetics 
of disordered vortex matter. Due to the complex na¬ 
ture of interacting vortex lines in the presence of disor¬ 
der, we need to carefully and systematically analyze the 
contributions of the different energetic contributions to 
fully understand this intriguingly complicated dynami¬ 
cal system, as well as assess the relevance of different 
initial configurations. Hence the implementation of tem¬ 
perature and magnetic field quenches in our numerical 
simulations constitutes a crucially important advance 
that should hopefully provide valuable insights for fu¬ 
ture experiments. Indeed, temperature or magnetic field 
quenches and the subsequent measurement of one-time 
observables and two-time correlation functions could be 
utilized to formulate experimental methods to character¬ 
ize material samples by means of the observed fluctua¬ 
tion spectrum. In future work, we plan to extend our 
study to include the effects of driving currents on the 
transient properties of driven flux line systems, where re¬ 
laxation towards non-equilibrium steady states will ensue 
in contrast to the presently analyzed relaxation towards 
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thermal equilibrium. A more detailed understanding of 
relaxation processes and the origin of deviation from scal¬ 
ing requires identifying the characteristic length scales in 
the system and their temporal growth, which we hope to 
pursue through the analysis of space-time observables. 
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